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Abstract 

We derive from first principles (without resorting to the replica trick) a density functional 
theory for fluids in quenched disordered matrices (QA-DFT). We show that the disorder-averaged 
free energy of the fluid is a functional of the average density profile of the fluid as well as the pair 
correlation of the fluid and matrix particles. For practical reasons it is preferable to use another 
functional: the disorder-averaged free energy plus the fluid-matrix interaction energy, which, for 
fixed fluid-matrix interaction potential, is a functional only of the average density profile of the 
fluid. When the matrix is created as a quenched configuration of another fluid, the functional can 
be regarded as depending on the density profile of the matrix fluid as well. In this situation, the 
replica-Ornstein-Zernike equations which do not contain the blocking parts of the correlations 
can be obtained as functional identities in this formalism, provided the second derivative of this 
functional is interpreted as the connected part of the direct correlation function. The blocking 
correlations are totally absent from QA-DFT, but nevertheless the thermodynamics can be 
entirely obtained from the functional. We apply the formalism to obtain the exact functional 
for an ideal fluid in an arbitrary matrix, and discuss possible approximations for non-ideal fluids. 
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I. INTRODUCTION 



The phase behaviour of fluids in quenched disordered matrices has been of prior interest 
in the last decade. The classical theoretical approach to these systems amounts to con- 
sider two different sets of state variables: the annealed variables (usually the position of 
the fluid particles), which are allowed to equilibrate, and the quenched variables (usually 
the position of the matrix particles), which have their values fixed. The reason for this 
distinction is that our system is not in thermal equilibrium with respect to the quenched 
variables, but it is in equilibrium with respect to the annealed variables for each fixed 
configuration of the quenched ones. Accordingly, two different statistical averages must 
be considered: the annealed average, which is the typical ensemble average of equilibrium 
systems, and the quenched average or average over disorder, which is performed over the 
quenched variables. For each configuration of the disorder, we can compute the equilib- 
rium thermodynamic magnitudes of the system by means of the corresponding annealed 
averages. These averages will, of course, depend on the configuration of the quenched 
variables. However, if the matrix is statistically homogeneous (its statistical features are 
similar everywhere) and the system is large, we expect little variation between annealed 
averages corresponding to different matrix configurations. Thus, quenched averages of 
the annealed averages are meaningful to characterize thermodynamic magnitudes of these 
systems. Because of this double average, the problem becomes intractable within the 
classical equilibrium statistical-mechanics tools and new theoretical methods are called 
for. 

Madden and Glandt j3| made an extension of the conventional diagrammatic treatment 
of liquid-state theory to obtain cluster expansions for the thermodynamics and structure 
of a fluid in a quenched matrix. They also derived a set of Ornstein-Zernike (OZ) equations 
relating the total and direct interparticle correlation functions, which can be solved with 
the appropriate closure relations p . Alternatively, Given and Stell Q| used the continuum 
version of the replica trick [4] to rederive this set of OZ equations and they noted that, 
although Madden and Glandt 's cluster expansions were correct, there were some missing 
terms in the OZ equations. The corrected set of OZ equations was called the replica 
Ornstein-Zernike (ROZ) equations and, for fluids with quenched-averaged density profile 
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Pi{x) and matrices with one-particle distribution po{x), is given by 

hoo{Xi,X2) =Coo{Xi,X2) + (CooPo ® hoo){Xi,X2), (1) 
hio{Xi,X2) =Cio{Xi,X2) + (CloPo ® hoo){Xi,X2) + (CcPi hio){Xi,X2), (2) 

/ioi(a3i,a;2) = coi(a;i, 332) + (cooPo ® /ioi)(iCi, 032) + (coipi (g) /ic)(a;i, a;2), (3) 
/iii(a;i,a;2) =cu{xi,X2) + (cioPo ® /ioi)(iCi, 0:2) + (ccPi ® 332) 

+ (cbPi ® /ic)(a3i,a;2), (4) 

hciXi,X2) =CciXi,X2) + (CcPl K)iXi,X2), (5) 

where (cp ® /i)(a;i,a;2) = / dccs c(a3i, a;3)p(cc3)/t(cc3, 0:2), the subscripts and 1 refer to 
the matrix and fluid, respectively, and 

hn{xi, X2) = hc{xi, X2) + hi,{xi, X2), (6) 
cn{xi, X2) = Cc{xi, X2) + Cb(a;i, 332), (7) 

where the subscripts c and b denote the connected and blocking parts, respectively, of the 
correlation functions. In terms of the replicated system, the blocking parts, /ib and Cb, 
are the zero-replica limit of the corresponding correlation functions between two different 
replicas of the fluid [3]. Clearly hio{xi,X2) = hoi{x2,Xi) and Cio(a;i,a32) = Coi{x2,Xi) 
and then it can be shown that Eqs. and (jHI) are equivalent, so that Eqs. (IH), 0, (jH) 
and © form an independent set. 

Rosinberg et al. 5| used the same replica trick to derive the thermodynamics of these 
quenched-annealed (QA) systems. There are two important results drawn from this work 
that concern the present paper: (i) the thermodynamics is completely determined by the 
connected parts of the correlation functions, and (ii) the connected and blocking parts of 
hii can be written without any reference to replicas as 



pi{xi)pi{x2)hc{xi, X2) = pn{xi,X2) - p{xi\{qi})p{x2\{qi}), 



Pi{xi)pi{x2)K{xi,X2) = p{xi\{qi})p{x2\{qi}) - pi{xi)pi{x2), (9) 
where denotes the quenched average, p{x\{qi}) is the equilibrium density profile of the 



fluid for a particular configuration {qi} of the disorder, pi(a3) = p{x\{qi}), and pii(a3i, X2) 
is the disorder-averaged pair correlation function of the fluid. 

The of Madde„ a„d G.andt Q, Give„ a„d SteU S and Ros^be., et al. S 

established the extension of the classical integral equation theory to fluids in quenched 



disordered matrices. Since then, this has been the main method to study QA systems 
and with its help much insight on the phase behaviour of these systems have been gained. 
But the rephca trick is closely linked to the ROZ equations and so it has the typical 
limitations of any integral equation theory: it is virtually impossible to apply the theory 
to non-uniform phases. In the case of fluids without disorder this problem was solved by 
density functional theories (DFTs), so it seems natural to ask for an extension of DFT to 
QA systems. 

There have been attempts to apply DFT to fluids in random media. For instance, 
Menon and Dasgupta have constructed a Ramakrishnan-Yussouff density functional, 
using the same replica trick employed in the derivation of the ROZ equations, to study the 
effect of pinning in the freezing of superconductor vortex lines. The same approach has 
been applied to study hard spheres in a quenched random gaussian potential More 
recently, Schmidt ^ has proposed a DFT for QA mixtures also based on the replica trick. 

In Schmidt's formalism, which we will refer to as replica-DFT (or simply rDFT), the 
matrix is described by the equilibrium free-energy density functional corresponding to the 
hamiltonian modelling the matrix particles, while the behaviour of the fluid is ruled by 
the quenched- averaged grand potential of the QA system ^tDft[Pi'i Po], which is written 
as a functional of the disorder-average density profile of the fluid, pi{x), and of the density 
profile of the matrix, po{q) (which enters as a parameter). The QA character becomes 
explicit in the minimization principle imposed over the quenched-average grand potential, 
which reads as 

= 0, (10) 

where po{q) is determined by the equation 

spoiq) 

FqIpo] being the equilibrium free-energy functional of the pure matrix and uo{q) = po — 
ipo{q), with pq the chemical potential of the matrix and fo{q) the external potential acting 
over the matrix particles. 

From i^rDFT[pi; Po], the free-energy functional can be defined as usual as 

^^rDFxbi; Po] = F'^^IPi] + i^rDFTbi; Po] - J dxUi{x) pi{x) , (12) 

where -F''^[pi] = kT f da[;pi(a:;)[ln Vipi(a;) — 1] is the ideal contribution (Vi being the 
thermal volume of the fluid particles) and Ui{x) = pi — ipi{x), pi being the chemical 



potential of the fluid and <fi{x) the external potential on the fluid particles. The excess 
contribution F^^pj,[pi; po] describes the interparticle interactions between fluid particles, 
and that between fluid and matrix particles. 

From a practical point of view, the functional -^rOFT [Pi ' Po] should be approximated. 
Schmidt's proposal for -FfDFT is based on fundamental measure theory P, [l^ This 
approximation has been applied to study the phase behaviour of colloid-polymer mixtures 
in bulk random matrices, of rods in quenched sphere matrices, of spheres in random fibre 
networks, and of soft-core fluids in soft-core matrices |l2]- Also, with the lattice version 
of fundamental measure theory it has been applied to study the freezing transition 
in a hard-core discrete fluid with different kinds of matrices ^|. Thus, as shown by its 
applications, rDFT is an important step forward in the study of QA systems. 

This notwithstanding, the theory has a number of weak points which should be pointed 
out. Although the replica trick is a widely-applied method of statistical physics, it makes 
a few assumptions which are difficult to justify concerning the analytic continuation of 
the grand potential as a function of the number of replicas, and the replica symmetry or 
its breaking. Hence an alternative derivation of DFT for QA systems would be desirable. 
Moreover, contrary to what happens in classical DFT, the formulation of rDFT makes it 
difficult to derive the set of OZ equations for QA systems from functional relations. As a 
matter of fact, at present it is not at all clear what the meaning of the second derivatives 
of -FrDFTlPi! Po] is- Thesc problems are the two main motivations of this paper. 

The remaining of the paper is organized as follows. In section 2, we propose a DFT 
for QA systems based on the convexity properties of the quenched-averaged grand po- 
tential. The derivation of the formalism resembles that of classical DFT and makes no 
use of the replica trick. We will Legendre-transform the grand potential to obtain the 
quenched-average "intrinsic" free-energy functional, J-'[pi, pio; Q], which depends on the 
quenched-averaged density profile of the fluid, pi{x), on the pair distribution function of 
the fluid and matrix particles, pio{x,q), and on the probability distribution of disorder 
Q. The dependence on piQ{x,q) can be eliminated to obtain a functional only of pi{x) 
(for fixed fiuid- matrix interaction), which will play the same role as the standard free- 
energy functional in classical DFT, and coincides with Schmidt's -FrOFT [pi ; Po] • We will 
then proceed with one of the most important contributions of this work: the derivation of 
a set of OZ equations where the direct correlation functionals are identified with second 
derivatives of this functional. Only equations (HHSI) and (0) can be derived within this 



DFT approach, but as we will discuss, these form a closed set of equations which involve 
all the structure information that is relevant to the thermodynamics of the system. Sec- 
tion 2 concludes showing how to derive the thermodynamics within this DFT approach. 
In section 3, we will make this formalism explicit in the case of an ideal fluid adsorbed in 
an arbitrary matrix. Conclusions and further discussions are gathered in section 4. 



II. QUENCHED- ANNEALED DENSITY FUNCTIONAL THEORY 

The system that we aim to describe consists of a fluid inside a porous matrix with which 
it interacts. The matrix is formed by a distribution of particles quenched at positions g,, 
i = 1, . . . , M. The fluid consists of particles, whose positions are denoted Xi, i = 1, . . . , N, 
and whose interactions are described by the hamiltonian (divided by kT) Ti^lxi, . . . , aj^v)- 
These particles are in equilibrium with a thermal bath at chemical potential fii and each 
of them undergoes the action of an external potential (fi{x). Besides, a fluid particle at 
position X interacts with a matrix particle at position q through the interaction potential 
ipio{x, q). To all purposes, the total external potential acting on a fluid particle at position 
X is 

M 

Veyit{x) = ^Piix) + (fwix, qi). 

i=l 

The grand partition function for this system will be {(3 = 1/kT) 



:[ui,Uio\{qi}] = 1 + ^ J , / dxi---dxN expi -H 

N=l ^ ' -J ^ 

N y M -| s 

+ Ui{Xi) + ^UiQ{Xi,qj) i, 

i=i L i=i -I J 



n[Xi, . . .,Xm) 



(13) 



which is a functional of ui{x) = fii — Lpi{x) and uiq{x, q) = —ipio{x, q), and also depends 
on the set {qi}- Accordingly, the grand potential will be 

n[ui,uw\{qi}] = -kTlnE[ui,Uio\{qi}]. (14) 

Now we need a model for the porous matrix. The simplest model is to assume that 
matrix particles are placed at random positions, according to a probability density Q|. 
Thus the grand potential is a random variable. The hypothesis we make now is that 
the grand potential per unit volume, in the thermodynamic limit, is a self-averaging 
random variable; therefore we can obtain its value in this limit by simply averaging over 



disorder (matrix particle positions). Hence the grand potential for the system in the 
thermodynamic limit is obtained as 

^l[ui,uio;Q] = ^[ui,Uio\{qi}], (15) 

where Q{{qi}) is the probability density of the matrix positions, and denotes a Q- 
average. This puts the quenched average into play. 



A. Concavity of the grand potential 

It is convenient to introduce the state functions 

N M 

Pn{x) = '^S{x - Xi), p^(q) = ^5(q- Qi). (16) 

i=l i=l 

In terms of them 

S [1*1,^10 1 {^i}] = 1 + ^ y^j^^ I dxi---dxN e^pl -nN{xi,...,XN) 

N=l 1 ■ I- 

+/3 (wi, Pjv) +/3(wio, PjvPm) I > (17) 

where 

{ui^pn) = j dxui{x)pN{x), {uio,pNp%) = J dxdquio{x,q)pN{x)p%{q). (18) 

With these definitions it is straightforward that 

5Q,[ui,uio; Q] 



p{x\{qi})p%{(i) =pio(a;,q), 



Suio{x, q) 

where p{x\{qi]) denotes the equilibrium density of the fluid for flxed positions of the 
matrix particles, and pi{x) is the quenched-averaged density profile of the fiuid. Likewise, 
Pio(a;, q) is the pair correlation function of the fiuid and matrix particles. 

On the other hand, VL[ui,uiq] Q\ is a concave functional of both, ui and uiq. This 
is easily proven by evaluating S[iii, iiio|{q'i}] on u 
(0 < A < 1). Since 

exp{ {uf^ ,Pn)}^ (exp{ {u^^^^ , Piv) }) (exp{ , Piv) }) 
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by Holder's inequality |3| we get 

S u[^\uio\{qi} < ^ 



(1) 

U\ ',Uio 



{q^} 



(0) 

U\ ',Uio 



{Q^} 



1~X 



from which 



U\ ',Uio 



{cl^} 



> XQ 



(1) 



{q,} +{i-x)n 



and averaging over disorder, 



> \Q 



u?,uio;Q\ +(l-A)fi \uf\uio;Q 



(0) 



{cl^} 



Clearly the same holds for uiq. Because of this, equations (fT^ define a one-to-one cor- 
respondence between the pair {pi,pio} and the pair {ui,uio} (i.e. the equations can be 
inverted) 



B. Free-energy functional and minimum principle 

Let us now introduce the Legendre transform of Q[ui,uio; Q] with respect to its two 
arguments 

^ [pi, Pio; Q] = ^[ui, uio; Q] + {ui, pi) + {uw, pw) , (20) 

where Ui{x) and UiQ{x,q) are the solution of Eqs. (11^1) for fixed pi{x) and pio(a;,q). 
Because of the properties of the Legendre transform [l^ 

(a) JF [pi, piQ] Q] is a convex functional of both pi{x) and pio{x, q); 

(b) the equilibrium pi{x) and pio{x, q) are the absolute minimum of the functional 



^ [pi, pio; Q]=0^ [pi, pio; Q] - {ui, pi) - {uw, pw) 



(21) 



for fixed iti{x) and UiQ{x,q), and therefore 
(c) they can be obtained by solving the equations 



^•^[pi,Pio; Q] 
5pi{x) 



ui{x), 



SJ^ [pi,Pio; Q] 

5pio{x,q) 



Uio{x,q) 



(22) 



As for the meaning of the functional JF, let us adopt a different point of view on the 
system: let us think of the porous matrix as an external potential acting on the fluid 
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particles. Then the "intrinsic" free-energy functional is obtained as 
T[p\ = VL[ui,Uio\{qi}] + j d.xp[x\{qi}) \^i{x) + ^^Uio{x, Qj) 

= VL[ui,Uia\{q.-^]+ / dx p[x\{qi})ui{x) + / da;dqrp(a;|{qi})p^^(q)Mio(a?, q) 



and is, of course, a functional of p{x\{qi\). If we now average this functional over disorder 
we obtain, making use of Eqs. (fT^ . (fT^ and (1201), 



^[p] = ^[pi,Pio;Q]. (23) 

This equation reveals the physical meaning of functional T[pi, pio; Q] as the intrinsic free 
energy of the fluid undergoing the presence of a porous matrix, averaged over disorder. 
But the fact that this functional depends on both pi{x) and pio{x,q) makes it rather 
inconvenient to use it as the basis for a QA-DFT (notice that this functional is the same 
for any fluid-matrix interaction, so it is far too general). 

Of course, one can assume uio{x, q) fixed and Legendre-transform only with respect 
to Ui{x) to obtain the alternative functional 

F[pi; Q] = Q[ui,uio; Q] + (Ml, Pi). (24) 

This is [for fixed Uio{x, q) and Q] a functional of pi{x) alone, and fulfils the Euler-Lagrange 
equation 

^ ^ (25) 
Comparing with and recalling that Uu){x,q) = —Lpio{x, q), 

F[p,; Q] = ^[pi, Pio; Q] + (<^io, Pio), (26) 

the intrinsic free energy of the fluid plus the interaction energy with the porous matrix. 
As we will show in Sec. Ill CI this one and not J^[pi,pio; Q] is the functional that plays a 
similar role in QA-DFT as the standard free-energy functional does in classical DFT, and 
in fact coincides with the functional -FrDFT[pi] derived from the replica formalism. 

C. Replica Ornstein-Zernike equations 

Let us work out the identity 

^ ' 5pr{x') J ^ 6m{y)6pi{x') J ^ 6u,{x)6uM 6pi{y)p,{x'y ^ ^ 
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This is one of the rephca Ornstein-Zernike (ROZ) equations, namely Eq. (0). To see it 
let us compute 

S'^n[ui,Uio; Q] 



-kT ' ' = pu{x,y\{q,i}) - pi{x\{qi})pi{y\{q,i}) - 5{x - y)pi{x\{qi}) 

dui[x)dui[y) 

= piix)pi{y)hc{x,y) - 6{x,y)pi{x), (28) 

where we have made use of Eq. (jSI). Introducing this expression into Eq. (P7|) we are 
immediately led to the identification 

dpi{x)pi{x') 

where F'^^[pi] Q] is the excess (over the ideal) part of the functional F[pi] Q]. 

In order to obtain Eqs. ((IHSl) we must set Q as the probability distribution of a grand- 
canonical ensemble at temperature Tq, chemical potential Pq and external potential ^Po{q)- 
Thus, if we define as usual Mo(q) = /^o ~ Vo{q)i the probability of finding the matrix 
configuration {M; qi, . . . , Qa/} is given by Q = {Puiqi, qj\/)}A/>o, where 

PM{qi,...,qM) = ^ ] 1 , .J, exp {-Hljiqi, qn) + (mo,Pm)} . (30) 
'='OFoJ Kq ^^J' - 

TiPj^{qi, . . . , qm) being the hamiltonian which models the interaction between matrix par- 
ticles (divided by /cTq) and So[uo] the grand partition function 

So[no] = 1 + ^ -^g— / dqi---dqMeicg{-nli{qi,...,qM) + {uo,p%)] (31) 

M=l ■ 

(Vo is the thermal volume of the matrix fluid). From classical DFT we know that for 
each external potential Mo(g) there exists a unique equilibrium density profile po(9) and 
the system can be described alternatively in terms of any of them. Therefore, if the 
hamiltonian Ti-M^qi, . . . , qn) and temperature Tq remain fixed, the dependence of the 
functional Q[ui] Q] and F[pi] Q] on the disorder is actually a dependence on either Uo{q) 
or po{q)- Hereafter, we will make explicit this dependence by writing this functional as 
F[pi; Pq] and the grand potential as Q[ui; Uq]. 
Now, Eq. ((H) is just the identity 

Spo{q') J 5uo{s) 5po{q') J 5uQ{q)Suo{s) 5po{s)po{q')' 

where fio[wo] = — ^7olnSo[po] is the grand potential of the matrix and -Fo[po] the corre- 
sponding free-energy functional. Finally, to obtain Eqs. and (jH)), we will notice that 

10 



the QA system can be described in terms of any of the following pairs of independent 
functions {ui{x) , uo{q)} or {pi{x) , po{q)} . Then, both ROZ equations can be identified, 
respectively, with the identities 
^ _ Sui{x) _ r^ 5ui{x)5pi{y) ^ r^5ui{x)Spois) 



Suoiq) J 5pi{y) 6uoiq) J 5po(s) Suo{q) 

6^F[p,;po] 6^n[u,;uo] 6'F[p,;po] ^^^oN 



dpi{x)6pi{y) 6ui{y)6uo{q) J 6pi{x)6po{s) 6uo{s)6uo{q)' 

(33) 



Spi{x) f ^^ Spijx) Suijy) r^Spi{x)Suo{s) 



Spoiq) J Sui{y) 6po{q) J 6uois) 6po{q) 



6ui{x)6ui{y) 6pi{y)6poiq) J 6ui{x)6uo{s) 6po{s)6po{q)' 

(34) 

where we have used Eqs. ()19p and (j^^ for the QA system, and their counterparts for 
the matrix. To complete the identification of these identities with the corresponding OZ 
equations, we have to take into account the expression ()28|1 as well as 

-^^/tiFtT = Pi{^)Po{q)h,o{x,q), (35) 
dui{x)duo{q) 

and to make the identification 

6pi{x)6po{q) 

At this point, it is important to notice that the set of OZ equations that we have 
obtained within the DFT approach is self-contained. This means that if we have the 
functionals -Fo[Po] s-nd F[pi] po], we can derive from them the direct correlation functionals 
Coo, Cc and cio, and using them as inputs in the OZ Eqs. ((T)), Q and (0) to obtain hoo, hio 
and he- Moreover, both Eqs. ^ and (jSJ can be solved independently and their solutions 
can be used to solve Eq. This situation is remarkably different from the one we 
find in the integral equation framework. In that case, although the OZ equation for the 
matrix [Eq. (0)] is independent of all the others, the remaining ones [Eqs. Q, (jlj) and 
(0)] form a coupled system. The reason for this difference is that in the case of integral 
equation theory, the direct correlation functions are also unknown and the ROZ equations 
must be complemented with closure relations. This additional equations are derived from 
exact relations between the interaction potentials between particles and the correlation 
functions, with one equation for each potential. Thus, in our case, we would have two 
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new equations which would involve {/iio,Cio} and {hu,Cii}, respectively. The absence of 
a closure relation for {he, Cc} is what keeps the set of Eqs. 0, (jH) and the two closure 
relations, and one of Eqs. (jHI) or ((Zj), coupled. 

We should remark that in the QA-DFT the blocking parts are absent. Nevertheless, 
contrary to what happens with integral equation theory, we are able to compute all 
the structure functions that are relevant to the thermodynamics without the blocking 
correlations. 



D. Thermodynamics 

We will finish this section showing how all the thermodynamics can be derived from 
the functional F[pi; pq]. The starting point will be the relation proved by Rosinberg et 
al. 



n[m\{q,}]=n[ui;uo] = -pV, (37) 
where p is the thermodynamic pressure and V the volume of the system (see Refs. j3l and 



17[ for a discussion about the definition of the thermodynamic pressure and the difference 
between this one and the mechanical pressure). Now, as the functional F[pi] po] is related 
to the quenched-averaged grand potential f2[Mi;Mo] [Eq. (j^ ] in the same way as the 
standard free-energy functional with the grand potential in classical DFT, and as this 
formal equivalence is also found in the relation ()25p between F[pi; po] and the chemical 
potential, we can conclude that all the thermodynamic relations we found in classical 
DFT remain formally identical in the QA-DFT, with Cc playing the role of the direct 
correlation because of Eq. (jSHl) • 



III. AN EXACT MODEL: IDEAL FLUID IN AN ARBITRARY MATRIX 

As the simplest example let us consider the only known example which can be exactly 
solved in this formalism: an ideal fluid in an arbitrary porous matrix. As a QA system, 
the matrix is taken to be a configuration of a grand-canonical ensemble of a certain fluid 
at temperature Tq, chemical potential po, and external potential <fo{q) [let us also define 
Uo{q) = pQ — ipo{q)]. The grand partition function, grand potential and free energy of this 
fluid will be denoted, respectively, Ho[uo], f^o[Mo] and -folPo], Po{q) being the corresponding 
equilibrium density profile. 
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For the ideal gas T^at = 0, so the grand partition function of the fluid becomes 

°° I ( r r ^ -I ^ ^ 

E[ui,Uio\{qi}] = 1 + < / da; exp Pui{x) + ^ puioix, q, 

M 



exp — I dx exp 



i=l 



(38) 



Thus, 



kT f 

n[ui,Uio\{qi}] = -— / da^ exp 



M 



Pui{x) + ^Puio{x, qi 



(39) 



and therefore 



Vi J ^oFoJ 
Vi 



(40) 



In this expressions Sq [^0(3?, ■)] stands for the grand partition function of the matrix fluid 
undergoing an external potential Uq{x, q) = UQ^q) + (To/T)uiq{x, q), x being the position 
of a fixed fluid particle, and AQo{x) = Qo[uq{x, ■)] — ^7o[wo]• 
From the first of Eqs. (fT^ it follows that 

, . 1 



PiW = — e 



f3ui(x)-l3oAno{x) 



(41) 



an interesting equation which tells us that the average equilibrium density profile of the 
fluid is given by the barometric law corrected with the probability of inserting a fluid 
particle in the matrix fluid at position x, namely e"'^"'^^''^^). 
From the second of Eqs. (fTn|) it follows that 

1 



(42) 



p,o(.,.) = ^e/^-^- L'f°/yiy"':^l 

Vi ^oKJ duo{x,q) 

Dividing this equation by Eq. (PT|) leads to 

Pioix,q) 6no[uo{x,-)] 
pi{x} 6uo{x,q) 

where po{x, q) is the equilibrium density profile of the matrix fluid corresponding to the 
external potential Uo{x, q) created by a fluid particle placed at x. 
Because of Eq. f)^T|) . Eq. simply becomes 



(43) 



^[ui,uio] Q] = -kT J dxpi{x) 



(44) 
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the equation of state of the ideal gas. On the other hand, ehminating Ui{x) from Eq. ()4H) . 



u,ix) = kT In (Vipi(x)) + —Ailoix), (45) 

so Eq. (pUj) becomes 

^[pi,Pio; Q] =F"'[pi] + ^/dx |pl(x)^]oK(^c, ■)] - Pl(a;)^^oN} 
+ J dx J dquioix,q)pio{x,q), 
where Mio(aJ, q) is the solution to Eq. and 

Fi'^lpi] = fcTy da;pi(x) |ln (Vipi(x)) - l} . (47) 



(46) 



Adding and subtracting 
T 

to Eq. (gni) and using Eq. (gSl) yields 



dxpi{x)!^j dquo{q)(^po{q) - po{x,q)^^ 



^[pi,Pio; Q] = F"^[pi] + j dxpi{x)S^j dquQ{q)(^po{q) - Po{x,q] 
+ dxpi{x) !^no[uo{x,-)] + j dquo{x,q)po{x,q] 



(45 



T 



dx pi{x) S^QqIuq] + j dquo{q)po{x] 



One can recognize in the brackets above the Legendre transforms of the grand potential 
of the matrix fluid; thus the final expression for the functional can be written as 



^[pi,Pio;Q] =F"'[p,] + ^J dajpi(a;){Fo[pio(a;,-)/pi(a;)]-Fo[po]} 
+ ^JdxJ dqr^^^|pi(a:;)po(g) -pio(a;,qr)|. 



(49) 



Notice that, apart from the standard ideal free-energy functional, there is a non-trivial 
term arising from the interaction between the fluid and the matrix. 

So far for the intrinsic free-energy functional. Now to obtain the functional F[pi; Q] — 
or, considering that we are describing a QA system, better F[pi; po] — we make a Legendre 
transformation of Q[ui,uio; Q] only w.r.t. ui {uio is assumed fixed) to obtain 



F[pi;po] =F'^[pi] + dxp^{x)[Fo[po{x, ■)] - Fo[po]} 

+ — / dxpi{x) / dq \pQ{q) ' - po{x,q)—^ \ 

To J J [ opo{q) opo{x,q) } 



(50) 
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A. The special case of an ideal matrix 

One particular case which has received some attention in the hterature 0, 0, Q] is 
the case in which the matrix is also ideal. The reason is that the ROZ equations for this 
system can be exactly solved (when uio is a hard-sphere potential), and, in spite of its 
simplicity, the blocking part of the direct correlation function is non-zero. 

If the matrix is a configuration of an ideal gas at temperature Tq, then 

FoIpo] = kTo j dqpo(q){ln (VoPo(q)) - l}- (51) 

Substituting this Fq in the expressions of the previous section one gets 

g/3oMo{q) 



Po(g) 



Vo ' 



„/3ouo{q)+/3«io(3;,q) 

Po(^, q) = ^ = po(g)e^"-(^''^); 



(52) 



therefore 



F[pi;po] = F"'[p,]-kT jdxj dqp,{x)po{q)Mx,q), (53) 
ho{x,q) = e/5-^o(-''?) - 1, (54) 



and 



g/3ui{a;) 



Pi(a;) = ^j^exp|y dqpoiq)fio{x,q)j . (55) 

A simple inspection of the exact functional ()53|) for an ideal fluid in an ideal matrix 
reveals that the non-ideal term is quadratic, so the only second derivative that is nonzero 
is Cio = /lo- Nevertheless, as mentioned above, this system has a non-trivial Cb 
which certainly cannot be derived from (|53p by any functional differentiation. In spite of 
this, the functional (j53|) contains all the equilibrium thermodynamics of the system. 



IV. DISCUSSION AND CONCLUSIONS 



We have made a first-principles derivation of a density functional formalism for 
fluids inside quenched disorder matrices without resorting to the replica trick. The 
main conclusion is that, for fixed interaction potential between the fluid particles, 
T-Cn{xi, . . . ,xn), and fixed distribution of the disorder, there exist a unique functional 
-^[PiiPioi Q] from which all the equilibrium structure information and thermodynamics 
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can be derived. Given the generalized external potential acting on the fluid particles, 
Ui{x) = fii — (fi{x), and the interaction potential between the fluid and matrix parti- 
cles, uio{x,q) = —Lpio{x,q), the disorder-average of the equilibrium density proflle of 
the fluid, p^^{x), and the fluid-matrix pair distribution p^^{x,q) can be derived from 
Eqs. ()22|) . Once we have Pi^{x) and p'^^{x,q), the average of the "intrinsic" free-energy 
of the system is given by J^[pi*^, Pio; Q], and the grand potential by 



Although for the ideal fluid in a quenched matrix we have been able to derive the 
explicit form of jF[pi,pio; Q], this is a formidable task for an arbitrary system. Note that 
this functional is valid for any interaction potential between the fluid and matrix particles 
and if we had it, then we would have solved a very general problem. Thus, it is more prac- 
tical to turn to a less general functional, F[pi; Q], which will be a functional only of pi{x) 
and whose functional form will depend on Uiq{x, q). Again, we have an Euler-Lagrange 
equation to obtain the equilibrium properties for a given generalized external potential 
Ui{x) [Eq. but now Flp^^; Q] is not just the average over disorder of the "intrinsic" 

free energy, but it also contains an additional contribution due to the quenched-average of 
the interaction energy between the fluid and matrix particles, — J dxdquiQ^x^ q)p'^{x, q), 
where p'^{x, q) can be obtained from the OZ Eq. (0). 

One of the most relevant contributions of this work is the identiflcation of the direct 
correlation functionals appearing in the ROZ equations with second functional derivatives 
of F[pi] Q] [Eqs. ()29|) and (jHUj) ]. It is worth mentioning that, in contrast to the case of 
classical DFT, the second derivative of F°''[pi; Q] with respect to pi{x) and pi{x') is not 
Cii(a3, x'; [pi; Q]) but only its connected part. Notwithstanding, the formalism is closed in 
the set of correlation functionals {(cqo, /ioo); (cio, ^lo)) (coi, ^oi)) (^c Cc)}, since the direct 
correlation functionals are obtained by simple functional differentiations of F[pi\ Q\ and 
the total correlation ones can be derived from the ROZ Eqs. (^jHl), and (0), which have 
been obtained as functional identities in the QA-DFT presented in this work. 

As we have discussed previously, the blocking parts /ib and Cb do not enter anywhere 
in the formalism. The fact that the QA-DFT does not contain these correlations and that 
the thermodynamics can be entirely derived from it implies that the blocking correlations 




PlO=Pio 



(56) 
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are not relevant for the thermodynamics. In this respect, we would like to mention 
that mode-coupling theory has been recently extended to QA systems in order to 
study the dynamics of confined glass-forming liquids, and the only equilibrium structural 
information needed to obtain the relaxing density fluctuations is the set {cc, Cio, Cq}. Thus 
even the liquid-glass transition can be determined if we know the functional F[pi; Q]. Also 
notice that in Refs. 0, Q| , where freezing is studied with a Ramakrishnan-Yussouff density 
functional, the direct correlation employed in its construction (which is derived with the 
replica trick) is Cc, not cn. 

Finally, as it happens in classical DFT, there are few systems for which F[pi; Q] can be 
obtained exactly (in this case only ideal fluids in arbitrary matrices, as far as we know). 
Therefore, this formalism should be complemented with approximations for F^^[pi; Q]. 

which make use of the 
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In this line are the works by Schmidt and collaborators 
constructing principle of fundamental measure theory [2, 
result for a OD cavity (a cavity which can hold at most either a fluid or a matrix particle) 
to approximate F[pi; Q]. Although the results obtained seem promising, we think that the 
extension of fundamental measure theory to QA systems involves subtleties concerning 
the correlations between fluid and matrix particles that are difficult to deal with, and 
further study is required. 

Another research line worth exploring is, in analogy to the development of classical 
DP approximations, to study the extension of those approximation based on the thermo- 
dynamics and structural information of the uniform fluid (usually obtained from integral 
equation theory) such as the weighted density or the effective liquid approximations j20|. 
This will be the subject of a forthcoming work. 
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